Quantum Persistence: A Random Walk Scenario 
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In this paper we extend the concept of persistence, weU defined for classical stochastic dynamics, 
to the context of quantum dynamics. We demonstrate the idea via quantum random walk and a 
successive measurement scheme, where persistence is defined as the time during which a given site 
remains unvisited by the walker. We also investigated the behavior of related quantities, e.g., the 
first-passage time and the succession probability (newly defined), etc. The study reveals power law 
scaling behavior of these quantities with new exponents. Comparable features of the classical and 
the quantum walks are discussed. 
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Persistence in classical dynamical systems is a topic 
which has been extensively studied during the recent 
years pj]. The persistence probability {Pci{t)) that the 
order parameter in a magnetic system has not changed 
sign till time t [11 and the persistence of unvisited sites 
in a diffusion problem Q are common examples which 
have received a lot of attention. The importance of per- 
sistence phenomena lies in the fact that the persistence 
probability in many systems shows an algebraic decay (in 
time) with an exponent not related to any other known 
static or dynamical exponents. 

The dynamics of a quantum system is expected to be 
different from the corresponding classical counter part, 
yet investigating persistence behavior in the quantum 
case has remained an interesting open problem till date. 
However, defining quantum persistence is ridden with the 
fundamental problem of measurement, since in order to 
ensure whether or not the system persisted in a given 
state (or in a chosen subspace X), one has to impose a 
continuous monitoring which would change the dynam- 
ics of a quantum system in some essential way. Hence in 
the quantum case, meaningful definition of persistence 
has to include the associated measurement scheme (i.e., 
how the evolution is disturbed by the measurement), and 
should essentially be dependent on that. The dynamical 
process we consider here is a discrete quantum random 
walk (QRW) 0-[i3- Classical random walk (CRW) on 
a line is a much studied topic where at every step 
one tosses a fair coin and takes a step, either to the left 
or right. The unitary implementation of QRW may be 
achieved through coupling an additional degree of free- 
dom (a quantum coin) with the walker. This coin degree 
of freedom is called the chirality, which takes values "left" 
and "right", analogous to Ising spin states ±1, and di- 
rects the motion of the particle. The state of the walker is 
expressed in the \x)\d) basis, where \x) is the position (in 
real space) eigenstate and \d) is the chirality eigen state 
(either "left" or "right, denoted by \L) and \R) respec- 
tively). There may be several ways of choosing the uni- 
tary operator causing the rotation of the chirality state. 



conventional choice effecting the rotation of the chiral- 
ity state is the Hadamard coin jsl-fioj unitary operator. 
(Most of the results are, however, believed to be coin 
independent.) The rotation is followed by a translation 
represented by the operator T: 
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The two component wave function il){x,t) describing the 
position of the particle is written as 
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and the occupation probability of site x at time t is given 
by 



f{x,t) = \i;L{x,t)\^ + \ijR{x,t)[' 



(3) 



normalization implying f{x,t) = 1. 

In the case of the classical random walker, two well 
studied quantities are persistence or the survival prob- 
ability Pci{x,t) defined as the probability that the site 
at X has not been visited till time t [1, 01 and the first 
passage time Fci{x,t) which is the probability that the 
walker has reached the site x for the first time at time t 
[6:]. The two quantities are related by 



F,i{x,t) = - 



dPci{x,t) 
dt ■ 



(4) 



At t >> 1, both the persistent probability and first pas- 
sage probabilities decay algebraically in time with expo- 
nents ad and /3ci which obey the relation 



(5) 



consistent with the equality in (|4]). We are primarily 
interested in the analogues of these two quantities in the 
QRW. 
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We now define our measurements schemes and the ob- 
servables to quantify the concept of quantum persistence 
in this case. In order to measure persistence in strict 
sense, one is left with no other choice than to monitor 
the system continually over time. One way of achieving 
this is to impose a direct time-continual projective mea- 
surement that determines at every moment whether or 
not the system persists within the subspace X in ques- 
tion. In this discrete time version of quantum walk, this 
amounts to carrying out a measurement after every time 
step of the unitary evolution following the scheme de- 
scribed below. The walk starts from some given site at 
t — 0, and a detector is placed at some other given site 
X, which detects the particle with probability unity if it 
reaches there. If the particle is detected at x, the evolu- 
tion is stopped (here, X is the entire lattice excluding x). 
Now the question asked for such a system (rather, for an 
ensemble of such systems), is what is the probability that 
the detector does not click till time t. This is the per- 
sistence probability P(x,t). It is equivalent to carrying 
out measurements at the site x after each step of unitary 
evolution of the ensemble and calculating the probability 
from the fraction of the surviving copies (for which x is 
yet unvisited) at each step. Within this setup of QRW 
on a line, placing the detector at x amounts to having 
a semi infinite walk (SIW) d UMl with an absorbing 
boundary at x and an open end in the other direction. 
Let us give a concrete illustration of the scheme with a 
detector placed at a; = 1. Suppose the walker starts at 
X = with left chirality. At time t = 1 in fifty per- 
cent cases it will be detected at x and the time evolution 
will be stopped. Persistence probability is therefore 1/2 
for X &t t = 1. The remaining fifty percent walks will 
evolve unitarily to the next step t ~ 2. At t = 3, the 
normalized probabilities aX x = —i^ x = —1 and x — 1 
are equal to 1/4, 1/2 and 1/4 respectively (and zero else- 
where). Hence now the detector detects the walker at 
X with probability 1/4, which means the 3/4 fraction of 
the population that was carried over to i = 2 would be 
carried over to the next time step at i = 3. This is 3/8 of 
the initial population (at t = Q). Hence the persistence 
probability at i = 3 will be 3/8 for x = 1. 

At each time step the ensemble is measured, and the 
amplitudes describe only to the surviving copies, and the 
probabilities are to be renormalized. Let the normal- 
ized occupation probability at x at time t be denoted by 
f{x,t). Thus f{x,t') denotes the fraction of the copies 
that survived the measurement at time t' — 1 (not the 
fraction of the initial population) which reaches x at time 
t' . The persistence probability is hence given by 

t 

Psiw{x,t)=l[{l-f{x,t'j). (6) 

t'=i 

It is to be mentioned here that by placing the detector 
at X, it is possible to find the occupation probabilities for 
all X and t (which are strongly dependent on x) but the 
persistent probability is obtained only for x = x. One 



may define a first passage time Fsiwi^jt) analogous to 
the classical random walk in this case as follows: 

t-i 

Fsiwix, i) = n (1 - 
t'=i 

= Psiw{x,t- l)f{x,t). 

It may be mentioned here that some related studies have 
been made earlier [ll|, [3 the problem of persis- 
tence measured in this way had been addressed with the 
boundary kept far from the starting point of the walker 

As mentioned before, definition of quantities in a quan- 
tum system depend heavily on the measurement scheme 
and we next pose similar interesting and well defined 
questions that brings out the more intrinsic character- 
istics of the dynamics somewhat directly, by monitoring 
what we call the succession probability S'(X, t) defined 
as follows. Let us consider a system allowed to evolve 
unitarily from a given initial state at ti = 0, up to a 
terminating time t' = At, when finally a measurement 
is done on it in order to determine whether or not it re- 
sides at a given state (or within a subspace) X and the 
evolution is stopped (e.g., in the QRW, one discards the 
walk). The entire process is repeated for increasing ter- 
minating times: t' = At, 2 At, 3 At, ...t. Now the question 
is asked, what is the probability that the system will be 
found within X in every measurement with t' < t. 

For a continuous-time evolution, this probability will 
be called the succession probability ^(X, t) in the limit 
At — > 0. For a discrete random walk. At will correspond 
to a single step. For example, in the context of QRW, 
one might choose to calculate the probability of a random 
walker not being found at some target site x in the suc- 
cessive measurements done at t' = l,2...t (starting from 
a given site). The subspace X in this case is constituted 
of all lattice points the walk may include, excluding x 
and we may then use the notation S{x,t) to denote the 
succession probability. It may be noted that in the clas- 
sical case, one need not restart the evolution, after each 
measurement, since the measurement would not disturb 
it. S{x,t) clearly differs from Pci{x,t) in general, since 
in case of S{x,t), in calculating the probability of find- 
ing the system within X at t', one takes contributions of 
all the paths running from the initial time up to the 
final time t' , including those which went out of X in the 
intermediate times. 

In the present one dimensional setting this amounts to 
allowing the system to evolve unitarily in either direction 
(infinite walk or IW [8]) for an interval t' , when the mea- 
surement is done and the walk is discarded. We choose, 
again, X to be the entire lattice except a given point 
X = X, where we would like to see whether the particle 
has reached or not. To determine S{x,t) for a given t, 
the termination time t' is varied as t' = l,2...t, and for 
each given t' one determines the occupation probability 
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f{x,t') and calculates S{x,t) as: 

t 

S{x,t)^l[{l-f{x,t')). (8) 

t'=i 

An analogue of the the first passage time may also be 
defined as 

t-i 

Fiwix, t)^Y[(l~ fix, t'))f{x, t) (9) 

= S{x,t-l)f{x,t). 

Experimentally this corresponds to simply knowing 
f{x,t') for t' < t. 

Some time dependent features in this type of infinite 
walk have been studied like the hitting time, recurrence 
time, Polya number etc. [l3 - [l7j . which involve the first 
passage time. However, in these studies, the spatial de- 
pendence has not been considered. For example, quan- 
tities like first passage time specifically at the origin has 
been estimated. 

It is important to note here that the quantities S and 
Psiw given by equations dS]) and dH) (as also Fjw and 
Fsiw) are identical in form: the difference being / ap- 
pearing in the infinite walk in place of / in the semi 
infinite walk. Thus / and / essentially make these quan- 
tities different. As an example, we have shown f{x,t) 
and f{x,t) as functions of 2; at a fixed time t in Fig[T] 
For the semi infinite walk, there is a detector placed at 
X = 10. To emphasize the difference, we have generated 
a walk biased towards the right, and the unbounded walk 
shows it clearly. On the other hand, in the semi infinite 
walk, the walker is not allowed beyond a; = 10 and con- 
sequently is driven towards the left. Obviously, even if a 
walk with symmetric boundary condition is initiated, the 
presence of the detector will convert it to a asymmetric 
walk. 

In the calculation, a quantum random walk is ini- 
tiahzed at the origin with -0^(0,0) = ao,ipji{0,0) = 
bo; a§ -I- 60 = 1. (All other Vl and tpn taken equal to 
zero). ipL{x,t) and ip}i{x,t) are recursively evaluated for 
all X and t. In the bounded (semi infinite) walk, contribu- 
tions from the walks going through x arc ignored. Unless 
otherwise specified, we have taken oq — 1/V2, bo — i/V^ 
which would result in a symmetric walk for the un- 
bounded (infinite) walk case. 

The results for SIW are essentially numerical; the per- 
sistence probabilities here saturate in time. This satura- 
tion behavior apparently originates from the simultane- 
ous effect of drifting of the quantum walker away from 
the origin and the presence of the boundary at x. These 
observations are in agreement with jl^ and consistent 
with other results involving recurrence time etc. (ol, [l3|. 
The first passage times, on the other hand, decay al- 
gebraically with t. As already mentioned, in [l^], the 
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FIG. 1: (Color online) Comparison of the probabilities f{x,t) 
for the semi infinite walk (SIW) and f{x,t) for the infinite 
walk (IW) at time t = 100: in the IW probabilities extend 
to both sides, in the SIW (with a detector placed at a; = 10), 
the particle is not found beyond x = 10. 




t/ Ixl 



FIG. 2: (Color online) The data collapse for the residual per- 
sistence probability in the quantum random walk shown for 
different values of a; — x. The straight line in the log-log 
plot has the slope = -2.0. Inset shows the variation of the 
saturation value Pq with x. 

persistence probability for large x was found to vary as 

Psiw{x,t) = Po +const(t/|a;|)-"^™' (10) 

where Pq is the saturation value and asiw — 2. We ver- 
ify this result with the observation that Pq has a weak 
dependence on x and observe that the numerical value 
of asiw approaches the value 2 asymptotically. The nu- 
merically estimated values of Po(x) are found to vary as 
(a — 6exp(— ex)) with c = 0.30 ±0.02, shown in the inset 
of FigO Using these values of Po(x), we show that a 
data collapse is obtained when the residual persistence 
probability Psiw{x.,t) — Po{x) is plotted against t/\x\. 
The first passage time Fsiw{x, t) behaves as 

Fsiw[x,t)^{t/\x\)-^^'^' /\xl (11) 

with Psiw — 3.0. Results for the collapsed data of per- 
sistence and first passage times are shown in Figs [5] and 
1 
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FIG. 3: (Color online) The data collapse for the first passage 
probability in the quantum random walk shown for different 
values of a; = a;. Fitted straight line has slope — -3.0. 
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FIG. 4: (Color online) The data collapse for the succession 
probability and scaled first passage probability defined for a 
infinite quantum random walk shown for different values of 
X ^ 10 (red), 20 (green), 50 (blue), 100 (magenta). 

For the unbounded or infinite walk, the calculations 
can also be done using the analytical forms available in 
i 

, / ,N ! + (-!)''+* f dk . cosfc , -Huj.t+kT 

^ 7 VI + cos^ k 

(12) 

2 J 2n + cos^ k 

(13) 

(which are obtained for a initial state with left chiral- 
ity, i.e., ao = l,bo = 0) and evaluate f(x,t) directly by 
numerical integration. 

Power law decay for both S{x,t) and Fiw{x,t) are 
observed: 

S{x,t) (X {t/\x\-l)-'^"^ (14) 

for t/\x\ > 1, and 

Fiwix,t)^{t/\x\^iy^'^/\x\ (15) 

for t/\x\ » 1 with aiw — 0.31 and Piw — 1-31- Data 
collapse for S{x, t) and Fjwix, t) from the numerical evo- 
lution of the infinite walk are shown in Fig. |4l These re- 
sults are obtained with = I/a/Z, = i/V^ which cor- 
respond to a symmetric walk. Results obtained from the 
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FIG. 5: Typical variation of Fiw{x,t) against x for different 
values of t = 10, 20, 50, 100, 200, 500. The curves extend over 
larger values of a; as t is increased. The peaks are approxi- 
mately fitted to 0.81a;;;°-|^. 

numerical integration of Eqs. [12] and [13] (corresponding 
to ao = 1 , &o = giving an asymmetric walk) show the 
same scaling behaviour as S{x,t) and Fiwix,t). Thus 
the exponents are independent of the initial conditions 
as expected. 

As discussed in [l^ [13] , the quantum walk on a line 
is recurrent, i.e. it returns to the origin with certainty 
and the same applies to visiting any other lattice point. 
Hence, the asymptotic succession probability is zero, 
which is in agreement with the power law decay of S{x, t) 
found in the present paper. 

While the exponents aiw ipisiw) and j3iw {l^siw) 
are different from the classical ad and /3c!, they enjoy 
a relationship identical to eq ([S]). We consider some 
other quantities related to the function Fiw{x,t) which 
one can compare with their classical counterparts. Plot- 
ting Fiwix^t) against x or we notice that it has 
an oscillatory behavior. These oscillations which die 
down for large values oi t/x as is apparent from Fig. 
[4] can be traced to the oscillatory behavior of f{x,t) 
for a QRW observed earlier |8l-[l0l|. From Figs [S] and 
[B] we observe that Fiw{x,t) actually attains a maxi- 
mum value Fiw,max{x,t) at values of = Xmax (or 
t — tmax) for fixed values of t (or x). We notice that 
Fiw,max{xjnax,t) oc x;;f^^ whcrc S ~ 0.59. Keeping x 
fixed, Fjw, max max 

) versus t^ax shows the same kind 
of dependence, i.e., Fiw,max{x,tmax) oc t"^^.. That the 
scalings with t^ax and Xmax turn out to be identical is 
not surprising as x scales as Mn a QRW. It is not possible 
to obtain this scaling form directly from eq. since 
F{x,t) attains a maximum value when t/\x\ is close to 
unity where the fitted scaling form is not exactly valid. 
In fact, eq. (|15p does not give any maximum value at all. 

Another dynamic quantity called hitting time has been 
estimated earlier for the QRW, in which an absorber is 
assumed to be located at a specific vertex of a hypercube 
within which the walk is conceived [13, [13 ■ The average 
hitting time is by definition the average time to reach that 
particular vertex for the first time. One can evaluate the 
average hitting time Th{x) in the infinite walk scheme 
using Th{x) — J2o tFiw{x,t) where t is allowed to vary 
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FIG. 6: Typical variation of Fiw{x,t) against t for different 
values of a; = 10, 20, 50, 100, 200, 500. The curves extend over 
larger values of f as a; is increased. The peaks are approxi- 
mately fitted to t^a^^ . 

from to T: 

Th ~ f tFjwix, t)dt 
Jo 

The numerical data (not shovifn) gives a fairly good agree- 
ment with this scaling. The above equation shows that 
Th blows up for T — > (X) in agreement with some earlier 
results using other coins [H, [3] ■ 

We thus observe that a number of quantities related to 
the dynamics of a quantum random walker follow power 
law behavior with time. Of these, the persistence proba- 
bility Psiw{x^t), which is obtained from a semi infinite 
walk, is drastically different from its classical analogue 
as it approaches a constant value in a power law fashion 
with an exponent which is quite different from the clas- 
sical value 1/2. The first passage time also has a power 
law decay with a new exponent. The numerical data also 
indicate that the two quantities obey a simple relation as 
in the classical case (eq (H))). 

A different quantum measure which we call the suc- 
cession probability has been proposed and calculated in 
the present work and a corresponding first passage prob- 
ability defined. These measures can be obtained from 
an infinite walk. The persistence probability and succes- 
sion probabilities are similar in form but the results are 
highly different with the succession probability exhibit- 
ing a power law decay (no saturation) with yet another 



new exponent. However, the form of the probabilities in 
eqs Uni), ([TT|) and (Ull), and the values of the expo- 
nents indicate the validity of eq (jU in the quantum case 
as well. 

In a classical random walk, {x^) cx V'^^ with 7c/ = 1 and 
in one dimension this scaling governs all other dynamic 
behavior including persistence. Thus all other exponents 
like ad and are essentially dependent on 7, e.g., ad = 
7c;/2 and Pd = f7c/- A quantum walker propagates 
much faster; here {x^) oc f'' with 7g — 2. Thus the 
dimensionless factor x/t appears in the scaling argument 
of the dynamic quantities. The exponents asiw — 2 and 
Psiw — 3 appear to be simply related to 7^; asiw — Iq 
and fisiw = §7g showing that here too the persistence 
phenomena is governed by the scaling (x^) oc t'^'' only. 

For the infinite walk case, the exponents aiw and j3iw 
are apparently not simply related to 7g. We have esti- 
mated some additional quantities involving the first pas- 
sage time. The maximum values of the classical probabil- 
ity Fd, behaves as l/tmax (for x constant) or l/a;^^^^^ (for 
t constant) showing that the obtained exponents are sim- 
ple multiples of ^d = 1. On the other hand, the behavior 
of Fiw,max appears to depend on the value of aiw and 
not 7g as it varies with tmax or Xmax with an exponent 
5 which is very close to 2aiw numerically. 

The average hitting time for a CRW is found to vary 
as T^"!-/"^. In the infinite QRW, this variation is given by 
rp2-i3iw _ Pqj. ^Yie classical case, 2 — f3d = "fcif^ but since 
no such relation exists for the quantum case, the hitting 
time scaling is therefore not dictated by 7^ but by Pjw 
(or ajw) only. 

Lastly, it is true that the probabilities f{x,t) are quite 
different from f{x,t) making the persistence and succes- 
sion probabilities distinct, however, the feature that the 
quantum walker walks away from the origin (in contrast 
to a classical walker) is present in both. This makes the 
persistence probability and the succession probabilities 
quite large in magnitude compared to classical persis- 
tence probabilities. For the quantum persistence proba- 
bility, Psiw 1 the additional constraint of the presence of 
the boundary makes it saturate. 
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